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Q '• Abstract 

\ We study the influence of errors and decoherence on both the performance of Shor's factoring 

> . 
O ' algorithm and Grover's search algorithm, and on the amount of interference in these algorithms 






using a recently proposed interference measure. We consider systematic unitary errors, random 

unitary errors, and decoherence processes. We show that unitary errors which destroy the interfer- 

Qh| ence destroy the efficiency of the algorithm, too. However, unitary errors may also create useless 
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Ch ' additional interference. In such a case the total amount of interference can increase, while the 
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efficiency of the quantum computation decreases. For decoherence due to phase flip errors, inter- 
ference is destroyed for small error probabilities, and converted into destructive interference for 



^ . error probabilities approaching one, leading to success probabilities which can even drop below the 

m ■ 

classical value. Our results show that in general interference is necessary in order for a quantum 



algorithm to outperform classical computation, but large amounts of interference are not sufficient 
and can even lead to destructive interference with worse than classical success rates. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Yz 



I. INTRODUCTION 

Quantum algorithms differ from classical stochastic algorithms by the facts that they 
have access to entangled quantum states and that they can make use of interference effects 
between different computational paths [l9|. These effects can be exploited for spectacular 
results. Shor's algorithm factors large integers in a time which is polynomial in the number 
of digits [4], and Grover's search algorithm finds an item in an unstructured database of size 
A^ with only ~ yN queries [5| . Many other quantum algorithms have building blocks similar 



^ 



% 18|. But more than twenty years 
91] it is still not clear what exactly is at 



to those developed in those two seminal papers (e.j 
after the discovery of the first quantum algorithm 
the origin of the speed up of quantum algorithms compared to their classical counterparts. 
Large amounts of entanglement must necessarily be generated in a quantum algorithm that 
offers an exponential speed-up over classical computation [lO|], and tremendous effort has 
been spent to develop methods to detect and quantify entanglement in a given quantum state 
12] for recent reviews). However, the creation of large amounts of entanglement 
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does certainly not suffice for getting an efficient quantum algorithm, and it remains to be 
elucidated what are both necessary and sufficient requirements. 

While there seems to be general agreement that interference plays an important role in 
quantum algorithms [l|, [j, [Sj, surprisingly, it has remained almost unexplored in compu- 
tational complexity theory. Recently we introduced a measure of interference in order to 
quantify the amount of interference present in a given quantum algorithm (or, more gen- 
erally, in any quantum mechanical process in a finite dimensional Hilbert space) [l3|. It 
turned out that both Grover's and Shor's algorithms use an exponential amount of inter- 
ference when the entire algorithm is considered. Indeed, many useful quantum algorithms 
start off with superposing coherently all computational basis states at least in one register, 
which is a process that makes use of massive interference (the number of i-bits, a logarithmic 



unit of interference, equals the number of qubits of the register 13|). Both algorithms differ 
substantially, however, in their exploitation of interference in the subsequent non-generic 
part: Shor's algorithm uses exponentially large interference also in the remaining part of the 
algorithm due to a quantum Fourier transform (QFT), whereas the remainder of Grover's al- 
gorithm succeeds with the surprisingly small amount of roughly three i-bits, asymptotically 
independent of the number of qubits. 



Recently it was shown that the QFT itself on a wide variety of input states (with efficient 
classical description) can be efficiently simulated on a classical computer, as the amount 
of entanglement remains logarithmically bounded [ij, [l5|, [l6| . As the QFT taken by itself 
creates exponential interference, it follows that an exponential amount of interference alone 
does not prevent an efficient classical simulation. This is in fact obvious already from the 
simple (if practically useless) quantum algorithm which consists of applying a Hadamard 
gate once on each qubit and then measuring all qubits. By definition, this algorithm uses 
exponential interference (T = 2" — 1 for n qubits). When applied to an arbitrary com- 
putational basis state, one gets any output between and 2" — 1 with equal probability, 
p = 1/2". According to Jozsa and Linden's result [lO(] this algorithm cannot provide any 
speed-up over its classical counterpart (as it creates zero entanglement), and indeed, it can 
evidently be efficiently simulated with a simple stochastic algorithm that spits out a random 
number between and 2" — 1 with equal probability, which can be done by choosing each bit 
randomly and independently equal to or 1 with probability 1/2. We note that the same 
phenomenon exists also for entanglement: the state (|000...000) + |lll...lll))/v2 has a lot 
of entanglement, but the corresponding probability distribution can be efficiently simulated 
classically with two registers. 

Thus the precise nature of the relationship between interference and the power of quan- 
tum computation is not yet fully understood. Surprisingly, there are tasks in quantum 
information processing which do not require interference in order to give better than clas- 



sical performance, as was shown in [17| for quantum state transfer through spin chains. 
In order to shed light on this question, in this paper we study both Grover's and Shor's 
algorithms in presence of various errors. We analyze to what extent the interference in a 
quantum algorithm changes when the algorithm is subjected to errors, and to what extent 
these changes reflect a degradation of the performance of the algorithm. We will investigate 
this question for systematic and random, unitary or non-unitary errors, and look at the 
"potentially available interference" as well as the "actually used interference", where the 
former means the interference in the entire algorithm, the latter the interference in the part 
of the algorithm after the application of the initial Hadamard gates [l3| . 



II. GROVER'S AND SHOR'S ALGORITHMS AND THE INTERFERENCE MEA- 
SURE 

As we will use Grover's and Shor's algorithms throughout the paper we first review 
shortly their main components. Grover's algorithm Uq [51] finds a single marked item a 
in an unstructured database of A^ items in 0{-\/N) quantum operations, to be compared 
with 0{N) operations for the classical algorithm. The algorithm starts on a system of 
n qubits (Hilbert space of dimension A^ = 2") with the Walsh-Hadamard transform W, 
which transforms the computational basis state |0...0) into a uniform superposition of 
the basis states A^~^^^ X^x^o 1^)- Then the algorithm iterates k times the same operator 
U = WR2WR1, with an optimal value k = [7r/(4^)] (where [. . .] means the integer part) 
and sin^^ = 1/N (isl, i.e. Ug = {W R2W Ri)''W . The oracle Ri multiplies the amplitude 
of the marked item a with a factor (—1), and keeps the other amplitudes unchanged. The 
operator R2 multiplies the amplitude of the state |0. . .0) with a factor (—1), keeping the 
others unchanged. 

Shor's algorithm ^ allows the explicit decomposition of a large integer number R into 
prime factors in a polynomial number of operations. The algorithm starts by applying 2L 
Hadamard gates to a register of size 2L where L = [logg -R] + 1 , in order to create an equal 
superposition of all computational basis states N'^/"^ X]t=o 1^) "where A^ = 2'^^ . Then the 
values of the function f{x) = a^ (mod R) , where a is a randomly chosen integer with < a < 
R, are built on a second register of size L to yield the state N'^/"^ St=o k)l/(^))- The last 
quantum operation consists in a Quantum Fourier Transform (QFT) on the first register only 
which allows one to find the period of the function /, from which a factor of R can be found 
with sufficiently high probability. In the numerical simulations, we didn't take into account 
the workspace qubits which are necessary to perform the modular exponentiation, but are 
not used elsewhere. It is a reasonable simplification in our case, since during this phase of 
modular exponentiation interference is not modified, as the whole process is effectively a 



permutation of states in Hilbert space [13|] . In order to study numerically the effect of errors 
for different system sizes, we performed simulations for n = 12 qubits, which corresponds, 
respectively, to factorization of i? = 15 (with a = 7), and also for n = 9 and n = 6. The 
cases n = 9 and n = 6 correspond to order-finding for R = 7 (with a = 3) and R = 3 (with 



a = 2) and do not exactly correspond to an actual factorization, although the algorithmic 
operations are the same, and a period is found at the end. In the case n = 9 the period 
found does not divide the dimension of the Hilbert space, so the final wave function is not 
any more a superposition of equally spaced 5-peaks, but is composed of broader peaks. This 
enables to reach the more usual regime of the Shor algorithm, where in general the period 
is not a power of two (although it does not happen for R = 15). In one case the result was 
different enough to warrant the display of the corresponding curve for a different value of a 
{a = 6) for which the period is a power of two (see end of Subsection III.C). 

The interference measure for a propagator P of a density matrix p {p[j = J2ki PijMPki) 
derived in [l3| is defined as 

^{P) = 2_^ \Pii,kl\ " Z-^ \Pii,kk\ ' ("'-) 

where Pij^ki are the matrix elements of the propagator in the computational basis {\k)} 
(fc = 0, . . . 2" — 1), and pki = {k\p\l). While the interference measure is certainly not unique, 
it quantifies the two basic properties of interference: the coherence of the propagation, and 
the "equipartition" of the output states, i.e. to what extent the computational basis states 
are fanned out during propagation. Indeed, the second term in ([1]) can be understood as 
a sum over matrix elements of a classical stochastic map (the map which propagates the 
diagonal matrix elements of the density matrix, thus the probabilities in the computational 
basis) . This term is subtracted from the more general first term, where the elements Pu^m of 
the propagator are responsible for the propagation of the coherences in the density matrix 
and their contribution to the final probabilities. Therefore, if all coherences get destroyed 
during propagation (i.e. the map is purely classical), interference vanishes. The squares in 
eq.([T]) are important, as they allow to measure the equipartition. The number of i-bits is 
defined as rij = log2(X(P)). One Hadamard gate provides one i-bit of interference jl3| . 
In quantum information theory the propagation of mixed states is generally formulated 



within the operator sum formalism 191]: A set of operators {Ei} acts on p according to 
p' = ^i EipEJ = Pp, where the Kraus operators Ei obey ^^^ ElE^ = 1 for trace-preserving 
operations. The interference measure then becomes ([1]) 
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In the case of unitary propagation presented by a matrix U, the interference measure reduces 
to 

J(P(f/)) = N-J2\Uik\'. (3) 

i,k 

This form makes it obvious that the interference is bounded by < I{P{U)) < N — 1. The 
interference measure is invariant under permutations of the computational basis states. 

III. PERTURBED QUANTUM ALGORITHMS 

In the following we examine the amount of interference in perturbed versions of Grover's 
algorithm and Shor's algorithm. We will distinguish between "potentially available" and 



"actually used" interference 13j. These names are motivated by the fact that both algo- 
rithms start from the single computational basis state |0), such that only the first column 
of the unitary matrix U, which represents the algorithms in the computational basis, deter- 
mines the outcomes and success probabilities. The interference measure, however, counts 
the interference for all possible input states, i.e. the contributions from all columns in U — 
thus the "potentially available" interference is in general much larger than what is needed 
in the algorithms. The actually used interference on the other hand is the interference in 
the remainder of an algorithm after all the initial Hadamard gates have been applied. At 
that point a coherent superposition of all computational basis states has been built up (in 
the case of Shor's algorithm: a coherent superposition of all computational basis states of 
the first register), and therefore all the interference measured by I{P{U)) is actually used. 
Another motivation to look at these two different measures is the fact that the latter focuses 
on the "non-generic" part of the algorithm. 

A. Systematic unitary errors 

We start by replacing each Hadamard gate with a perturbed gate, parametrized with an 

angle 9, as 

(cos 9 sin 9 \ 
sin 9 — cos 9 I 
The unperturbed Hadamard gate corresponds to ^ = 7r/4, while the cases 9 = and 9 = 7i/2 
replace the Hadamard gates by the Pauli matrices a^ and a^, respectively, which create no 
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FIG. 1: (Color online) Potentially available interference in the Grover algorithm with systematic 
unitary errors in the Hadamard gates, parametrized by the angle 0, eq. ([!]) (left); success probability 
5* of the algorithm (middle); and success probability as function of interference (right). Black circles 
mean n = 4, red squares n = 5, green diamonds n = 6, blue triangles up n = 7, gray triangles left 
n = 8. All curves are averaged over all values of the searched item a. 

interference. Thus the replacement of Hadamard gates by (j4]) amounts to destroying the 
interference produced in the course of the algorithm in a controllable fashion. This allows 
us to compare the loss of interference with the efficiency of the algorithm in a systematic 
way. We measure this efficiency through the success probability S. For Grover's algorithm 
the natural defintion of S is the probability to find the searched state |q;), S = tr(p/|a)(a|), 
where the density matrix pf describes the final state at the end of the computation, which 
may be a mixed state if decoherence strikes during the calculation (see section UlI CI) . For 
Shor's algorithm there are in general many "good" final states which allow to compute the 
period of the function /, and it is therefore more appropriate to define 5* through the loss 
of probability on these "good" final states compared to the unperturbed algorithm. Thus, 
if X^i'^/'iK) is ^^^ final state of the unperturbed algorithm and Yli'^T^l^)^ ^^ define S for 
Shor's algorithm as 



5 = i-Eli^^i'-i^ri'|/2. 



(5) 



Figure [H shows the dependence of the potentially available interference and of the success 
probability 5* of Grover's algorithm on 6, as well as the success probability as function of the 
interference. All curves are averages over all values of a, a = 0, . . . , 2" — 1. Both interference 
and success probability peak at 9 = n/A. For a small number of qubits, S{9) shows some 
additional modulation in the wings of the n curve, which are pushed further and further 
out for increasing n. The broad maxima of X(^) lead to steep increases of S{X) close to the 
maximum possible value for the interference X = 2" — 1. At^ = Oor0 = 7r/2, the interference 
vanishes, as in that case the algorithm degenerates to a combination of permutations and 
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FIG. 2: (Color online) Same as Fig. [H but for actually used interference. 
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FIG. 3: (Color online) Potentially available interference in the Shor algorithm with systematic 
unitary errors in the Hadamard gates, parametrized by the angle 6, eq. ^ (left); success probability 
S of the algorithm (middle); and success probability as function of interference (right). Black circles 
mean n = 6 (/(x) = 2^ (mod 3) ), red squares n = 9 (/(x) = 3^ (mod 7) ), green diamonds n = 12 
(/(x) = 7^ (mod 15) ). Inset on the left is a close-up of the case n = 12 close to the maximum. 

phase shifts, so that no two computational states get superposed. Figure [1] shows that for 
this example an exponential amount of interference is necessary even in order to obtain a 
success probability of the order 1/2, and by squeezing out a small additional amount of 
interference, 5* can be boosted to its optimal value. 

The actually used interference gives a similar picture: For 6 = or 6 = 7t/2 the inter- 
ference vanishes, and interference reaches its maximum value T ^ 4 for 6 = n/A. As the 
success probability remains unchanged whether we calculate the interference for the entire 
algorithm or only after the initial Hadamard gates, we find again that the success probability 
increases with increasing interference (see Figj2l). 

Figure [3] shows the potentially available interference and the success probability for Shor's 
algorithm. Again, both interference and success probability peak at 6 = tt/A. The additional 
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FIG. 4: (Color online) Same as Fig. [3l but for actually used interference. 



modulation in the wings of the curve for S{9) for a small number of qubits is much less 
pronounced than for Grover's algorithm, but the main fact remains that the broad maximum 
of 2{6) corresponds to a sharp peak for S{9), leading to the same steep increase of S{X) 
close to the maximum possible value of the interference. However, the exact algorithm does 
not lead to the maximum possible amount of interference. Close to ^ = 7r/4, the interference 
slightly increases while the success probability goes down, indicating that some interference 
which is "useless" in terms of the algorithm efficiency is generated. When 6 increases and 
the interference is reduced by a large amount, the algorithm has a low success probability. 
Even though the success probability globally goes down with an increasing number of qubits, 
it is not the case for all 6 values. This can be attributed to the fact that the three curves 
on the figure do not exactly describe the same problem for different numbers of qubits, but 
are instances of order-finding for different values of the number R. Thus when different 
values of the number of qubits are used, the precise problem investigated depends on the 
number theoretical properties of the integers chosen, which can be different. On the contrary, 
Grover's algorithm run on different numbers of qubits is essentially the same problem run 
on a computer of different size. 

In (iSj, it was pointed out that Grover's and Shor's algorithm use a very different amount 
of actually used interference. Indeed, for Grover's algorithm it remains bounded for all values 
of the number of qubits n , while for Shor's algorithm it grows exponentially with n. This 
may be related to the fact that Shor's algorithm is exponentially faster than all known 
classical algorithms, while Grover's is only quadratically faster. In Figure H] the actually 
used interference is plotted for Shor's algorithm, showing that for 6 = 7r/4 it reaches its 



fi- ^ - o OOJ00O00000OO( 




'■ [ I 

ill} 



FIG. 5: (Color online) Potentially available interference in the Grover algorithm with random uni- 
tary errors in the Hadamard gates, parametrized by the interval e, eq.(|l]) (left); success probability 
S of the algorithm (middle); and success probability as function of interference (right) for n = 4 
to n = 7. Same symbols as in Fig. [TJ 

maximal value which grows exponentially with the number of qubits. In this case, any 
decrease in the interference corresponds to a decrease of the success probability. 

B. Random unitary errors 



Let us now consider what happens if we replace each Hadamard gate with a gate given by 
eq.(jl]), where each 6 is chosen randomly, uniformly and independently from all other gates 
in an interval it /A — e/2, 7r/4 + e/2. 

Figure [5] shows the interference and success probability of Grover's algorithm as function 
of e. All curves are averaged over all possible values of a as well as over rir random realizations 
of the algorithm {rir = 1000 for n = 4, rij. = 100 for n = 5, ... ,7). Again, the maximum 
amount of interference is obtained for the unperturbed algorithm, e = 0, but the maximum 
is not very prominent. It is followed by a shallow minimum close to e = 7r/4, which gets 
shifted to smaller values for increasing n. Altogether, the interference is little affected by 
the random unitary errors. This can be understood from the fact that the unitary matrix U 
representing the algorithm is already almost full in the unperturbed algorithm 13|, with the 
exception of the first column, which propagates the initial state |0), and presents a strong 
peak on the searched item. Randomly replacing the Hadamard gates by H{6) increases the 
equipartition in the first column, as is witnessed by the decay of 5*, but reduces on average 
the equipartition in the other columns, leading to a slight overall decrease of interference 
close to e = 0. This means again that a very large amount of interference is necessary in 
order to get even a modest performance of the algorithm, and a very steep increase of the 
success probability occurs when interference is boosted to its maximum value. 
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FIG. 6: (Color online) Actually used interference in the Grover algorithm with random unitary 
errors in the Hadamard gates, parametrized by the interval e, eq.dl]) (left), and success probability 
S as function of interference (right). Same symbols as in Fig.[TJ The success probability as function 
of e is the same as in Fig. [5J 



The situation for the actually used interference is quite different, as shows Fig. [61 For 
e = 0, we have an interference X ~ 4 at the end of the algorithm (after the first diffusion 
gate it reaches its maximum value of X ~ 8 — 24/ N and then oscillates and decays with each 
subsequent diffusion gate to the final value X ~ 4 [l3|). Thus, the unperturbed algorithm 
leads to remarkably low equipartition in the entire matrix U, a highly unlikely situation 
for any random matrix. Indeed it was shown in 20(] that a random unitary N x N matrix 
drawn from the circular unitary ensemble (CUE) gives, with almost certainty, an interference 
X ~ A^ — 2. Thus, it is not surprising that with growing e, X rapidly increases to a value 
I r^ N . As the success probability decreases with e, this leads to the counter-intuitive 
situation that the success probability decays with increasing actually used interference. 

Figures [TfHl display the effect of random unitary errors on Shor's algorithm with the 
number of random realizations n,. = 5000 (n = 6), n^ = 1000 {n = 9), and Ur = 100 
(n = 12). Besides changing the Hadamard gates, random phases with the same distribution 
were added to the two-qubit gates in the quantum Fourier transform. This way in both 
algorithms, all Fourier transforms and Walsh- Hadamard transforms are randomized in a 
comparable way. Figure [7] shows that the potentially available interference oscillates slowly 
as a function of e, on a scale which seems independent of the number of qubits and also 
larger than for Grover's algorithm. The situation is similar to the one in Fig. [6] (although 
Fig. [6] deals with actually used interference), since interference increases for small values of 
e and reaches a maximal value around e = 0.6 — 0.7, while the success probability decreases. 

The situation is different in the case of actually used interference, shown in Fig. [HI Indeed, 
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FIG. 7: (Color online) Potentially available interference in Shor's algorithm with random unitary 
errors in the Hadamard gates, parametrized by the interval e, eq.dH) (left); success probability of 
the algorithm (middle); and success probability as function of interference (right). Symbols as in 
Fig. [3l Inset shows a close-up of the curve for n = 12. 
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FIG. 8: (Color online) Actually used interference in Shor's algorithm with random unitary errors 
in the Hadamard gates, parametrized by the interval e, eq.([l]) (left), and success probability as 
function of interference (right). Same symbols as in Fig. [3l Inset shows a close-up of the curve for 
n = 12. The success probability as function of e is the same as in Fig. [71 

interference starts from its maximum possible value and decreases with increasing e at the 
same time as the success probability decreases. In the same way as for potentially available 
interference, the variation of interference is relatively small compared to the case of sys- 
tematic errors, which were explicitly designed to destroy interference. Nevertheless, Fig. [8] 
shows that contrary to the case of Grover's algorithm, interference and success probability 
decrease in a correlated way. 
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C. Decoherence 

We finally consider a class of errors which create true decoherence. We distinguish be- 
tween phase flips and bit flips, and consider a (somewhat artificial) situation, where the 
errors occur only during the first Walsh-Hadamard transformation, i.e. the sequence of 
Hadamard gates on all qubits at the beginning of the algorithm in the case of Grover's 
algorithm, and on all qubits of the first register of length 2L in the case of Shor's algorithm. 
We will assume that nj out of n qubits are affected by errors, and study interference and 
success probability as function of n^, nj = 1, . . . , ra. Note that if all Hadamard gates in the 
entire algorithm were prone to error, one would need to calculate 2^^'^+^^" Kraus operators for 
Grover's algorithm (see secjlT]), each of which is a 2" x 2" matrix, which makes the numerical 
calculation rapidly too costly. The former number is reduced to a more bearable 2^^ in our 
case. For Shor's algorithm, the number of Hadamard gates depends on the implementation 
of the modular exponentiation and the calculation of the function /, but grows exponentially 
with the number of qubits as well, if all qubits can be affected by the decoherence process. 
Contrary to the calculations for unitary errors, in the simulation of Grover's algorithm we 
restricted ourselves to a fixed value of the searched item a, but checked for a few different 
values of a that the results are insensitive to the value of a. 

A Hadamard gate prone to errors is followed with probability p by a bit-fiip (or by a 
phase flip — we consider only one type of error at a time), and we calculate again both 
the potentially available and actually used interference. In the latter case, only the Pauli 
operators y^a^ and ^Jp(yx which represent the phase flip and bit flip errors, respectively, 
with probability p on a given qubit are included in the Kraus operators, but not the initial 
Hadamard gates themselves. 

For Grover's algorithm. Fig. [9] shows the result in the case of bit flip errors, and Fig. [10] 
for phase flip errors. For both types of error, the interference has maximal value for p = 
oi p = 1, which corresponds to completely coherent propagation, and minimal value for 
p = 0.5. The minimal value decreases rapidly with the number of qubits prone to error. The 
potentially available interference reaches zero in the case of bit flip errors on all n qubits, 
whereas for phase flip errors a finite value remains. The actually used interference shows the 
opposite behavior. It has zero minimal value for phase flip errors on all n qubits, whereas it 
remains finite for bit flip errors even with probability 0.5. Phase flip errors rapidly destroy 



13 



FIG. 9: (Color online) Potentially available interference in the Grover algorithm with decoher- 
ence through bit-flips during the first Walsh-Hadamard transformation, as function of the bit-flip 
probability p after each Hadamard gate (left). Same but for actually used interference (middle). 
Success probability as function of p (right). All curves are for n = 4, a = 2; nj = 1 black circles, 
Uf = 2 red squares, nj = 3 green diamonds, Uf = A blue triangles. 




FIG. 10: (Color online) Same as Fig. [9l but for phase flip errors. 



the operability of the algorithm. The success probability decreases linearly with p for nj = 1 
to reach a value close to zero for p = 1 , and more and more rapidly for increasing nj . In fact , 
for n = 4, S{p = 1) ~ 0.0025, independent of nj, which is even smaller than the classical 
value 1/16=0.0625. The algorithm is completely coherent in this case, and the large amount 
of interference is used in a destructive way, subtracting probability from the searched item. 
Remarkably, bit-flip errors do not affect the success probability at all, such that S{p) 
remains constant at the optimal value, independently of the number of qubits affected. 
The behavior is easily understood, as in fact the bit flip errors leave the state obtained after 
applying the Hadamard gates invariant (and in particular: pure). The interference goes down 
to zero, nevertheless, as it measures coherence using superpositions of all computational 
states 13|, and not pureness of the flnal state. We have therefore the peculiar situation 



where in spite of decoherence processes one particular state remains pure (the perfectly 
equipartitioned superposition of all computational basis states), and since it is that state 
which is used in the algorithm, the success probability remains unaffected. On the other 
hand, the interference measure was constructed to measure coherence by the sensitivity of 
flnal probabilities to relative initial phases between the computational basis input states, 
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FIG. 11: (Color online) Potentially available interference in the Shor algorithm with decoherence 
through bit-flips during the first Walsh-Hadamard transformation, as function of the bit-flip prob- 
ability p after each Hadamard gate, for n = 12 (left), n = 9 (center), n = 6 (right). The symbols 
are: nj = 1 black circles, Uf = 2 red squares, nj = 3 green diamonds, nj = 4 blue triangles up, 
Uf = 5 cyan triangles down, nj = 6 brown stars, Uf = 7 gray x, n/ = 8 violet +. The correspond- 
ing success probability 5 is constant equal to p = 1 for all values of p (data not shown). Here and 
in the following figures all quantities are averaged over all possible choices of the nj qubits in the 
first register. 

and it correctly picks up that the phase coherence between all states got lost. Thus, in this 
particular situation, one can have a perfectly well working algorithm which uses, according 
to our measure, zero potentially available interference. We believe, however, that this case 
where the coherence of the propagation cannot be measured by the influence of relative 
phases but only through the purity of the final state is highly exceptional and should not 
exclude a proof that exponential speed-up needs exponential interference, if one restricts 
attention to unitary algorithms, or gives special attention to the exceptional single pure 
state mentioned. It is also important to note that the actually used interference remains 
finite at p = 0.5 and nj = n, and below the already small value T ~ 4 for the unperturbed 
algorithm and thus never goes to zero for bit-flip errors. 

Figures [TT1IT21 show the result of decoherence due to bit-flips in the Shor algorithm for 
n = 12, n = 9 and n = Q qubits. As for the Grover algorithm, bit-flips are performed after 
each of the initial Hadamard gates. However, as these Hadamard gates concern only one of 
the registers, the decoherence process affects only the first two-thirds of qubits. The curves 
show the effect of decoherence on a growing number of qubits, from nf = 1 to Uf = 2/3n, 
with data averaged over the choice of the nf affected qubits. The success probability S is 
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FIG. 12: (Color online) Actually used interference in the Shor algorithm with decoherence through 
bit-flips during the first Walsh-Hadamard transformation, as function of the bit-flip probability p 
for n = 12 (left), n = 9 (center), n = 6 (right). The corresponding success probability S is constant 
equal to S* = 1 for all values of p (data not shown); symbols as in Fig. [TTl 



computed as in eq.(l5]) where now iV'f^'^P is replaced by the probability for state i in the final 
mixed state. 

The success probability is constant equal to 1 for all values of p (data not shown). In 
contrast, both potentially available and actually used interference are strongly affected by 
the decoherence. Both quantities decrease from their maximum value at p = and p = lio 
a minimum at p = 0.5, the potentially available interference decreasing faster. However, the 
interference never goes down to zero in this setting, as can be seen in the insets of Figs llllllSl 
contrary to the case of the Grover algorithm above. 

In the Figs JTTKTSl only two-thirds of the qubits are affected by the decoherence. This may 
appear to be the main reason why the potentially available interference does not decrease to 
zero in presence of bit-flip decoherence, contrary to the case above with the Grover algorithm. 
In order to investigate in more details this question, we studied the interference produced 
when all n qubits are affected by bit-flip decoherence. The results are shown in Fig. [13] for 
n = 6. Although the two types of interference decrease faster than in Figs JTTl[T2l none of 
them reaches exactly zero over the whole interval of p values. Contrary to the case of FigsITTl- 
T2l the success probability is now strongly affected by the decoherence and is not preserved: 
the initial superposed state is not any more protected against this type of decoherence, since 
the second register is not supposed to be in a equal superposition state in the exact algorithm. 
The data displayed in Figs. fTTIITS] suggest that for Shor's algorithm, although it is possible to 
decrease the interference by a large amount while keeping the success probability constant, 
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FIG. 13: (Color online) Interference and success probability in the Shor algorithm with decoher- 
ence through bit-flips during the first Walsh-Hadamard transformation, as function of the bit-flip 
probability p after each Hadamard gate, with all n qubits flipped, for n = 6: potentially available 
interference (left), actually used interference (center), success probability (right). Inset on the left 
is a close-up close to the minimum; symbols as in Fig. [TTJ 




FIG. 14: (Color online) Potentially available interference in the Shor algorithm with decoherence 
through phase-flips during the first Walsh-Hadamard transformation, as function of the phase-flip 
probability p after each Hadamard gate, for n = 12 (left), n = 9 (center), n = 6 (right). Inset on 
the left is a close-up of the case n = 12 close to the minimum; symbols as in Fig. [TTJ 

it does not seem possible to perforin efficiently the computation with zero interference. 

Figures [I^lfl6l display data obtained for decoherence through phase-flips. As before, phase 
flips are introduced with probability p on qubits of the first register after application of the 
Hadamard gates. The data displayed in Figs JTUfT^ show that interference, both potentially 
available and actually used, decreases to a minimum at p = 0.5. The minimum is lower 
than in the case of bit-fiips, and reaches zero for actually used interference. In the case of 
potentially available interference, some residual interference is still present when all qubits of 
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FIG. 15: (Color online) Actually used interference in the Shor algorithm with decoherence through 
phase-flips during the first Walsh-Hadamard transformation, as function of the phase-flip proba- 
bility p after each Hadamard gate for n = 12 (left), n = 9 (center), n = 6 (right). Inset on the 
left is a close-up of the case n = 12 close to the minimum, showing that the value X = is indeed 
reached for p = 0.5. 



s 

0.8 * 

0.6 - >; 

0.4 - 

0.2 - 





,N.- 



0.2 0.4 0.6 
P 



*s-g-tf 6" 6 6-0 Q ffl a 6 




FIG. 16: (Color online) Success Probability in the Shor algorithm with decoherence through phase- 
flips during the first Walsh-Hadamard transformation, as function of the phase-flip probability p 
after each Hadamard gate for n = 12 (left), n = 9 a = 3 (center left), n = 9 a = 6 (center right), 
n = 6 (right); symbols as in Fig. [TTJ 

the first register are affected. Figure [16] shows that in contrast to the case of bit-flip errors, 
success probability is strongly affected for phase-flip errors. It decreases with the number 
of qubits affected and the value of p until the algorithm is totally destroyed. Comparison 
with the Figs JTHfTSl shows that the increase of interference between p = 0.5 and p = 1 
is not reflected in a similar increase in success probability, the interference produced in 
this case is useless and does not serve the algorithm. It is similar to the one found for 
random algorithms in [20|], where it was shown that random algorithms produce on average 
an interference close to the maximum value. The case n = 9, a = 3 is peculiar: in this 
particular instance (second figure in Fig. [T6l) . after an initial decrease for small values of 
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p, the success probability increases for larger values of p, although it never reaches values 
close to one. We think this is due to the fact that in this case the period is not a power of 
two, and therefore the final wavefunction is composed of broad peaks which for such small 
sizes have a significant projection on many basis states of the Hilbert space. A random- 
type wavefunction produced by the destroyed Shor algorithm therefore has a much larger 
projection on such state than on a state composed of sharp S peaks as in the two other 
values of n. To sustain this hypothesis, we computed the success probability for n = 9 and 
a = 6, where the period does divide the Hilbert space dimension, and the final wave function 
is composed of (5-peaks. In this case (last figure of Fig. [161) the success probability indeed 
goes to zero for large p values. 

IV. CONCLUSIONS 

In this paper we have investigated the success probability of quantum algorithms in re- 
lation with the interference produced by the algorithms. To this end we subjected Grover's 
search algorithm and Shor's factoring algorithm to different kinds of errors, namely sys- 
tematic unitary errors, random unitary errors, and decoherence due to bit flips or phase 
flips. The study of systematic unitary errors showed that in both algorithms the controlled 
destruction of interference goes hand in hand with the decay of the success probability. This 
reinforced the idea that interference is an important ingredient necessary for the functioning 
of these algorithms. 

The case of random unitary errors shows, however, that a large amount of interference is 
by no means sufficient for the success of a quantum algorithm, since in some cases the inter- 
ference increases with decreasing success probability. This effect is particularly pronounced 
for the actually used interference in the case of Grover's algorithm, where the interference 
increases from about two i-bits for the unperturbed algorithm to an amount of the order 
n i-bits, close to the maximum possible value, for sufficiently strong errors. The success 
probability may decrease even below the classical value corresponding to unbiased random 
guessing, meaning that the interference has become destructive. This can also be under- 
stood in the context of the recent result that a randomly chosen quantum algorithm leads 



with very high probability to an amount of interference close to the maximum value |20|, 



such that randomizing a given quantum algorithm with limited interference is very likely to 



19 



increase the interference, even if the algorithm itself is destroyed in the process. Thus, as to 
be expected, interference needs to be exploited in the proper way to be useful. 

The study of decoherence led to more complex results. Phase-flip errors destroy interfer- 
ence and in parallel decrease the success probability, in the same way as systematic errors. 
Interference decreases for small error probabilities, and reappears again as destructive in- 
terference for error probabilities approaching one, leading to success probabilities which can 
even drop below the classical value. In contrast, bit-flip errors performed after each initial 
Hadamard gate do not reduce the success probability of the algorithm, while affecting the in- 
terference produced. In the case of Grover's algorithm, the potentially available interference 
can even go all the way down to zero while the performance of the algorithm is unaffected. 
This surprising result is due to the symmetry of the equipartitioned state used as initial 
state in the Grover's algorithm, which is invariant under bit-fiips. It should be remarked, 
however, that the actually used interference does not go to zero for Grover's algorithm. As 
concerns Shor's algorithm, the bit flips destroy part but not all of the interference, both 
potentially available and actually used, while also keeping constant the success probability 
of the algorithm. These results show that in general it is possible to reduce substantially 
the interference produced while keeping the efficiency of the algorithm. Grover's algorithm 
seems to run correctly with zero potentially available interference, but not with zero actually 
used interference. In contrast, in none of our simulation Shor's algorithm was found to run 
efficiently without some interference left, potentially available or actually used. 
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